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Characterising Body Tissue 



Field of the Invention 

The present invention relates to methods for the characterisation of body 
tissue. More specifically, the invention is concerned with the characterisation 
of body tissue as normal (e.g. healthy) or abnormal (e.g. pathological). The 
invention has particular, although not necessarily exclusive, applicability to the 

r- - 

diagnosis and management of cancer, including breast cancer. 
Background 

In order to manage suspected or overt breast cancer, tissue is removed from 
the patient in the form of a biopsy specimen and subjected to expert analysis 
by a histopathologist. This information leads to the disease management 
program for that patient The analysis requires careful preparation of tissue 
samples that are then analysed by microscopy for prognostic parameters such 
as tumour size, type and grade. An important parameter in tissue 
classification is quantifying the constituent components present in the sample. 

> 

Interpretation of the histology requires expertise that can only be learnt over 
many years based on a qualitative analysis of the tissue sample, which is a 
process prone to intra observer variability. 

* 

* - 

Despite the relative value of histopathological analysis, there remains a 
degree of imprecision in predicting tumour behaviour in the individual case. 
Additional techniques have the potential to fine-tune tissue characterisation to 
a greater degree than that currently used and hence will improve the targeted 
management of patients. 

In existing research in this field, x-ray fluorescence (XRF) techniques have 
been used to study trace element composition of breast tissue and have 
shown that breast cancer is accompanied by changes in trace elements and 
such measurements could contribute to tissue grading 1 . It has also been 
shown that x-ray diffraction effects can operate as an effective means of 
distinguishing certain types of tissue 2,s . Furthermore, it has been shown that 
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such diffraction effects could be suitably analysed to demonstrate small 
differences in tissue components and that this analysis could lead to a 

A 

quantitative characterisation of tissues . 
Summary of the Invention 

5 It is a general, preferred aim of the present invention to develop quantitative 
analytical approaches that could add precision to the characterisation of 
tissues, in particular to distinguish between normal and diseased (e.g. 
pathological) tissue. 

A preferred aim is to add precision to the several subjective components of 
10 tissue analysis, most notably those variables 'scored' in breast tumour 
grading. 

In general terms, the invention provides methods for analysing and/or 
characterising body tissue in which results are obtained by considering a 
combination of two or more different types of measured tissue properties. 

15 In one aspect, the present invention provides a method analysing body tissue, 
the method comprising: 

obtaining data representing a first measured tissue property of a body 
tissue sample; 

obtaining data representing a second, different tissue property of the 

20 tissue sample; and 

using the data in combination to provide an analysis of the tissue 

sample. 

In another aspect, the invention provides a method for characterising body 
tissues as normal or abnormal the method comprising: 
25 obtaining data representing a first measured tissue property of a tissue 

sample; 

obtaining data representing a second, different tissue property of the 

tissue sample; and 

using the data in combination to provide a characterisation of the tissue 

30 sample as normal or abnormal. 
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In either aspect it is preferred that data representing a third measured tissue 
property is also used in combination with the other data in the analysis or 
characterisation of the tissue sample. 

It is particularly preferred that data representing four or more measured tissue 
5 properties is used in combination in the analysis or characterisation of the 
tissue sample. 

Suitable techniques that can be used to obtain the tissue property data 
include x-ray fluorescence (XRF), energy or angular dispersive x-ray 

> 

diffraction (EDXRD), Compton scatter densitometry, low angle x-ray scattering 
10 and the measurement of linear attenuation (transmission) coefficients. 

The tissues properties that are measured may include the composition of the 
tissue sample, for instance the presence, concentrations and/or proportions of 
specific elements or organic compounds. 

Preferably, in either of the aspects above, the data is used in combination to 
15 obtain the desired result by using the data as the input to a predefined 
calibration model that relates the combined data to one or more tissue 
characteristics, (e.g. normal or abnormal). 

In a further aspect, the invention provides a method for creating a tool for the 
analysis and/or characterisation of body tissue, the method comprising 
20 creating a calibration model that relates data representing two or more 

(preferably three or four or more) measurable tissue properties to one or more 
tissue characteristics. 

The calibration model is preferably produced by using sets of the measured 
data from tissue samples for which the characteristic(s) (e.g. normal / 
25 abnormal) to be determined by the model are already known. These data 
sets can be used to 'train 1 the model in a known manner. 

Other multivariate analysis techniques may be employed. 
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Brief Description of the Drawings 

Embodiments of the invention are described below by way of example with 
reference to the accompanying drawings, in which: 

Figure 1 is a schematic diagram of EDXRD experimental apparatus employed 
in the exemplary methods described below according to embodiments of the 
invention; 

Figure 2 is a series of graphs showing average XRF responses; 
Figure 3 shows EDXRD scatter profiles for normal and diseased tissue; 
Figure 4 shows PLS model predictions for the normal test samples;, 
Figure 5 shows PLS model predictions for the diseased test samples; 
Figure 6 shows predictions of tissue type for the normal test samples; and 
Figure 7 shows predictions of tissue type for the disease test samples; 
Description of Embodiments " 

The following describes a study that has been conducted in vitro using breast 
tissue samples. The principles exemplified below are, however, more widely 
applicable, for example to other tissues and in relation to data obtained from 
in vivo as well as in vitro measurements. The principles can usefully be 
employed in in vivo imaging applications. 

The exemplary embodiments described below employ x-ray fluorescence 
(XRF) and energy dispersive x-ray diffraction (EDXRD) techniques to reveal 
the tissue characteristics. The invention is not, however, limited to these two 
techniques and other techniques may be used in addition or as alternatives to 
XRD and EDXRD. Other techniques that might be used include Compton 
scatter densitometry, low angle x-ray scattering and linear attenuation 
(transmission) coefficients. 
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As discussed in more detail below, during the study the concentrations of K, 
Fe, Cu and Zn were measured in 77 breast tissue samples (38 classified as 
normal and 39 classified as diseased) using X-Ray Fluorescence (XRF) 
techniques (in other embodiments, concentrations of other elements or 
5 organic compounds might be measured). The coherent scattering profiles 
were also measured using Energy Dispersive X-Ray Diffraction (EDXRD), 
from which the proportions of adipose and fibrous tissue in the samples were 
estimated. 

The data from 30 normal samples and 30 diseased samples were used as a 
10 training set to construct two calibration models, one using a Partial Least 
Squares (PLS) regression and one using a Principal Component Analysis 
(PCA) for a Soft Independent Modelling of Class Analogy (SI MCA) technique. 
The data from the remaining samples, 8 normal and 9 diseased, were 
presented to each model and predictions were made of the tissue 
15 characteristics. 

Three data groups were tested: XRF, EDXRD and a combination of both. The 
XRF data alone proved to be most unreliable indicator of disease state with 
both types of analysis. The EDXRD data was an improvement, however with 
both methods of modelling, the ability to predict the tissue type most 
20 accurately was by using a combination of the data. 

1. Breast Tissue Samples 

The tissue samples measured were obtained from mastectomies, 

r 

lumpectomies and breast reduction surgery. In regard to the latter, a number 
of healthy breast tissue samples were obtained. The tissue obtained from 

25 mastectomies or lumpectomies, was generally taken from the site of a lesion, 
classified as invasive ductal carcinoma, and in some cases normal tissue was 
taken from areas distant to the tumours. In line with the available samples, 
investigations were made for 38 samples classified as normal and 39 samples 
classified as diseased. The weight of each of the specimens was of the order 

30 of 1 g. Most specimens were of thickness in the range of 2-3 mm. Following 
excision the samples were kept frozen at - 85° C, no processing or sample 
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preparation taking place between excision and measurement. For both the 
XRF and the EDXRD measurements the samples were allowed to thaw 
before being measured in room temperature. 

2. Experimental Procedure 

2. 1 XRF 

t 

The XRF studies were carried out making use of the European Synchrotron 
Radiation Facility (ESRF), working on the Bending Magnet beamline BM28 5 . 
Using a simple arrangement of incident synchrotron radiation, tuned to a 
photon-energy just above the K-absorption edge of interest, particularly low 
elemental detection limits are achievable (< 1ppm). The high intensities of 
XRF available allow for short measurement times, providing for a high sample 
throughput. For the synchrotron photon beam, the plane of polarisation is the 
same as that of the electron orbit. Thus for a 90° geometry between photons 
directed on to the sample and the normal to the detector (Si(Li), Gresham 
Scientific Instruments, Sirius model), the strong linear polarisation of the 
photon beam provides significant suppression of the scattered photon 
intensity (fluorescence being unaffected). Given that the detector has lateral 
extent, the remaining sample-dependent scattered photon (coherent and 
incoherent) intensity reaching the detector is therefore governed by the solid 
angle formed between the sample and detector crystal. In addition to 
providing improvement in the signal to background ratio, control of the 
scattered radiation intensity allows use of the scattered peak area as a 
normalisation factor. As tissue is a low Z material, the fact that the detection 
system cannot resolve the Compton component will not affect the results. 

Each element of interest (K, Fe, Cu and Zn) was identified by the photopeak 
associated with its K a fluorescence photon emission. In order to seek L 
maximum production of the K a photons of interest, the samples were 
irradiated by photons of energy 500 eV above the particular K absorption 
edge, being an arrangement which also allows for the resolution of the 
scattered incident peak and the fluorescence response. The exception to this 
method was for K, where the data were collected using the same incident 
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photons as that for Fe. In order to quantify the sample concentrations of the 
elements of interest, calibration curves were constructed for each element. 

The calibration standards were aqueous solutions of the elements, the water 
matrix of the calibration models matching the "wet" nature of the tissue 
specimens. The following ranges of concentrations were used for the 
calibrations, as indicative of those expected to be found in tissue: 

K: 100, 300 and 1000-4000 ppm in increments of 1000 ppm 

Fe: 3-30 ppm in increments of 3 ppm 

Cu: 1-10 ppm in increments of 1 ppm 

Zn: 2-25 ppm in increments of 2 ppm 

The calibration solutions were measured in petri dishes that were sealed with 
laboratory sealing film (LabSeal, Merck). The tissue specimens were placed 
on such petri dishes previously filled with purified water and sealed. The 
specimens were then covered with the same sealing film. The beam size on 

* 

the specimens was 3 mm x 0.5 mm. 

The spectra acquired from the standard solutions were analysed using the 
software PeakFit (PeakFit™ SPSS Inc, AISN Software lnc)„developed for 
spectroscopy. The spectra were smoothed using a procedure based on * 
deconvolution, leading to the removal of peak broadening effects caused by 
imperfect resolution of the measuring instruments. The spectra were 
subsequently fitted using a procedure based on the Levenburg-Marquardt * 

* ■ 

non-linear minimisation algorithm. The fitting process took into account a 
linear baseline resulting in an estimation of the net total of counts integrated 
over the width of the photopeak. In order to normalise the fluorescence 
response, the scattered photopeak area was also calculated. The ratio of 
fluorescence to scattered photon peak area was then used to derive the 
relationship between element fluorescence and its concentration. 
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The tissue samples were irradiated under the same conditions as those used 
for the standard solutions and spectra were collected for each of the elements 
of interest. The analysis of the spectra also followed the procedure described 
for the standard solutions. The least squares fit derived from the calibration 
data, which relates the ratio of fluorescence to scatter photon peak area and 
. the element concentration was used to quantify the levels of each element in 
each of the samples. It is acknowledged that only a small area of the sample 

k 

is irradiated but measurements indicate that the inhomogeneity has been 
found not to significantly alter the profound differences between healthy and 
cancerous tissues. No correction has been made for matrix effects has been 
made in this study. However, the interest here is in the comparison between 
the levels of healthy and cancerous tissue and as any errors are systematic 
the comparison is not compromised. 

2.2 EDXRD 

» 

The EDXRD scatter profile of each of the samples was measured using a 
technique that utilises the scatter of a^polyenergetic photon beam at a fixed 
scatter angle. This technique has been used for a number of biomedical 
applications, notably that of estimating bone mineral density 6,7 and more 
recently for breast tissue analysis 2 . For an overview of applications of X-ray 
diffraction analysis in crystalline and amorphous body tissue see as an 
instance Royle et al 8 . 

The experimental set up is shown schematically in figure 1 . A tungsten target 
x-ray tube operating at 70 kV and 15 mA was used, the intrinsic filtration being 
1 mm of beryllium. The incident beam is tightly collimated via a slit cut in a 
dural slab forming a rectangular cross section of dimensions 1 mm x 2 mm. A 
similar collimation arrangement was set up at a scatter angle of 6° leading to 
a scattering volume in which the thickness of the sample was enclosed. The 
samples were mounted in standard 35 mm slide frames, sealed on either side 
with film 4 microns thick (Ultralene from Glen Spectra Reference Materials). 
The translator enabled the samples to be moved through the incident beam, in 
this instance a distance of 3 mm, producing an approximate irradiated volume 
of 12 mm 3 . The scattered photons were detected using an HpGe detector 
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(EG&G Ortec). The output pulses were analysed using a multi channel 
analyser (92X Spectrum Master, EG&G Ortec). The measurement time for 
each sample was 2400 seconds. 

3. Results 

3. 1 XRF 

Figure 2 shows the average of all XRF spectral response for the four elements 
of interest for normal and diseased tissue. Quantitative values of the 
elemental concentrations were obtained from the ratio of the XRF response 
peak to the scattered peak (not shown in figures) via the calibration line. 

Table 1 summarises these results for all tissue samples, showing the range of 
concentrations within a sample group, the mean concentrations and the 
standard deviation. For details about the measured levels of the elements 
and the comparison between healthy and cancerous tissue see work by 
Geraki et al 1 . 



Units In 
PPM ' 


K 


Fe 


Cu . 


Zn 


Normal 


Diseased . 


Normal 


Diseased 


Normal 


Diseased 


Normal 


Diseased 


Range 


57-514 


114-1120 

to 


< D.L.- 
36.64 


< D.L.- 
48.98 


< D.L.- 
1.53 


< D.L- 
2.09 


< D.L.- 
7.51 


< D.L.- 
20.19 


Mean 


163 


512 


8.82 


17.62 


0.27 


0.88 


2.19 


6.96 


SD 


141 


336 


8.74 


14.09 


0.22 


0.52 


1.37 


5.31 



Table 1 : Summary of the XRF results for the elements K, Fe, Cu and Zn. 
3.2 EDXRD 

Figure 3 shows the averaged diffraction spectra for all normal and all diseased 
tissue samples. The difference in the composition of the two types of 
specimens is evident. 

The characteristic peak from the adipose tissue can be seen at a momentum 
transfer value of 1 .1 nnrf 1 and the characteristic peak from fibrous tissue is at 
approximately 1 .6 nm" 1 . These peaks were fitted using the same technique as 
for the XRF spectra. The evaluated photopeak areas reflecting the presence 
of adipose and fibrous tissue were then corrected for factors such as the 
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shape of the x-ray tube spectrum and the difference in the attenuation and 
scatter properties of the two types of tissue. The corrected relative intensities 
of the two scatter peaks reflect the relative amounts of the two materials in 
each specimen. 

These results show a strong dependence of the relative amounts of the two 
materials on the type of specimen. The healthy specimens were 
predominately made of fat (76 ± 9 %) while the tumour specimens were 
mainly composed of fibrous tissue (85 ± 4 %). 

4. Modelling 

The above data were divided into two groups i.e. 30 normal and 30 diseased 
sample data were used as a training set to produce a calibration model. The 
data from the remaining samples from each group (8 normal and 9 diseased) 
were then used as input for the model and the tissue type was predicted. 

In the calibration process we use the empirical data (i.e. elemental 
composition and EDXRD data) and prior knowledge (i.e. quantities known to 
represent that data) to predict unknown quantitative information from future 
measurements. In this simple application we are using a multivariate 
approach i.e. we use many variable measurements X{x(1), x(2)....x(n)} to 
quantify the target variables Y{y(1), y(2)....y(m)}. In this case the X variables 
are the data measured above and the Y variables are the categories of tissue 
i.e. normal or diseased. 

In nrr\&r to nre^rWrti the* l^ter frnm th^ former \a/a n^Rrl tn ^Qtim^tA how Y mtatA?; 

to Y and an example of this is a regression model of the form y=Xb+c. The 
unknown parameters b and c are estimated from the calibration data which 
can then be used to predict future values of y from measured X. 

Another method is to use a classification procedure where models are created 
that represent a particular classification of variable. This is carried out using 
principal component analysis on the data belonging to a particular 
classification. Input data are then analysed and compared with the models 
and a fit to each classification is established. 
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There are several statistical calibration methods available for creating a 
model 9 and this work uses a partial least squares (PLS) regression and 
principal component analysis (PCA), which is used for Soft Independent 
Modelling of Class Analogy (SIMCA) classification. The Unscrambler 10 is a 
software package that was used to perform multivariate analysis in this study. 

4. 1 Partial Least Squares (PLS) regression 

Three models were created using the PLS approach. The first used the XRF 
data alone, the second used the EDXRD data alone and the third used the 
combined data set. The appropriate data from the test samples (8 normal, 9 
diseased) were then presented to the model and predictions of tissue type 
made. 

Figures 4 and 5 show the predictions for the normal samples for all three 
models and the predictions for the diseased samples respectively. 

4.2 Classification models 

For each of the data groups (XRF, EDXRD and combined data) principal 
component analysis models were constructed for the healthy and the 
cancerous samples. The appropriate data from the test samples were 

presented to the models and a score was obtained indicating how close each 

,— 

sample came to each of the models. 

ij 

Figure 6 and 7 show the predictions for each sample in the normal category 
and diseased category respectively. 

5. Discussion 

It can be seen that the XRF model predictions using the PLS approach are the 
most unreliable with the EDXRD and the combined data being similar. 
However the uncertainty in the predictions becomes significantly smaller when 
both data sets are combined. 
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If an acceptable prediction parameter was chosen at 70% certainty and 
above, the number of reliable and false predictions can be found. Table 2 
summarises the data for the PLS model 





XRF 


EDXRD 


Combined 




Normal 


Diseased 


Normal 


Diseased 


Normal 


Diseased 


Mean 

Prediction 

(%) 


72.8 


69.8 


82.2 


98.5 


82.6 


99.2 


Mean 

Uncertainty 

i%) 


11.2 


32.9 


20.1 


21.6 


12.2 


14.5 


True +ve 


3 


2 


3 


6 


5 


7 


False -ve 


0 


1 


0 


0 


0 


0 


Undecided 


5 


6 


5 


3 


3 


2 



5 Table 2 : The mean predictions, uncertainty and number of true, false and 
undecided predictions for each data group and tissue type using PLS. 

Similarly, when using the classification technique, the XRF is the most 
unreliable with an improvement being shown using the EDXRD data. It 
should be noted that for the normal samples predictions were higher in the 
10 normal classification whereas for the diseased samples the wrong 

classification was made in several instances. The use of the combined data 
shows a marked improvement in prediction particularly when examining the 
diseased samples. 

As above, choosing a 70% probability cut off limit, the number of accurate 
15 predictions can be found. Table 3 summarises the data. 





XRF 


EDXRD 


Combined 




Normal 


Diseased 


Normal 


Diseased 


Normal 


Diseased 


Mean 

Prediction (%) 


66.8 . 


53.2 


67 


59 


77.8 


77.2 


True +ve 


3 


0 


4 


3 


6 


8 


False +ve 


0 


1 


0 


1 


0 


0 


Undecided 


5 


8 


4 


5 


2 


1 



Table 3 : The mean predictions, true and false positives for each data group 
and tissue type using the classification technique. 

The relative inefficiency of the XRF data (compared to the combined XRF and 
20 EDXRD results) in successfully and accurately predicting the type of the test 
specimens is due to the wide spread of concentrations that characterise the 
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groups of samples, evident by the large associated standard deviations (table 
1). 

As illustrated by the study described above, embodiments of the present 
invention can provide improved characterisation of tissue types using a 
combination of data and an appropriate model. A classification technique has 
been shown to be particularly successful. 

Embodiments of the invention have been described above byway of example. 
It will be appreciated that various modifications to that which has been 
specifically described can be made without departing from the invention. For 
instance, the study described above to exemplify the invention involved the 
use of only two types of tissue characterising properties. Other embodiments 
of the invention may use more than two types of tissue characterising 
properties or alternative characterising properties. Creating a model using 
samples that are characterised using a variety of useful parameters may 
develop useful histopathology tools. Provided the different data groups can 
represent all the parameters one wishes to quantify, the multivariate approach 
is a promising method for accurate characterisation of samples. 
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Figure 1 

Schematic diagram of the EDXRD experimental apparatus. The primary and 
scatter slit collimators are 1.0 x 2.0 mm and the scatter angle is 6°. 
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Figure 2 

Shows the average XRF response for K, Fe.Cu and Zn. Note also the 

presence of Ar in air and Ca in tissue. 
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Figure 3 



EDXRD scatter profiles for normal and diseased tissue, (average of 

samples) 
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Fiaure 4 



The PLS model predictions for the normal test samples for each data type i 
XRF, EDXRD and the combination of both. The error bars represent the 

uncertainty in the predictive ability. 
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Figure 5 

The PLS model predictions for the diseased test samples for each data type 
i.e. XRF, EDXRD and the combination of both. The error bars represent the 

uncertainty in the predictive ability. 
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Figure 6 

Predictions of tissue type for the normal test samples using the classification 

technique for the three data groups 
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Figure 7 

Predictions of tissue type for the disease test samples using the classification 

technique for the three data groups 
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